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Abstract. Orbiting x-ray and gamma ray instruments are snbject to large background count rates due to local 
particle fluxes in the space environment. The ability of an instrument to make calibrated measurements of the flux 
from a source of interest is highly dependent on accurately determining the background level. We present here a 
method of calculating the energy dependent background flux for any point in the complete data set recorded by 
the Burst and Transient Source Experiment (BATSE) in its nine year mission. The BATSE Mass Model (BAMM) 
uses a Monte Carlo mass modelling approach to produce a data base of the gamma ray background which is then 
Hltered to simulate the background count rate with a 2.048 second time resolution. This method is able to reduce 
the variations in the background flux by a factor of 8 - 10, effectively ‘flat-fielding’ the detector response. With 
flat-fielded BATSE data it should be possible to use the Earth occultation technique to produce a hard x-ray all 
sky survey to the 1-2 mCrab sensitivity limit. BAMM is also capable of estimating the contribution to the spectra 
measured from gamma ray sources due to the reprocessing of source photons in inactive material surrounding a 
gamma ray detector. Possible applications of this aspect of the model in the area of Gamma Ray Burst spectral 
analysis are discussed. 

Key words. Gamma rays: observations - Space vehicles: instruments - Methods: data analysis - Background 
simulation 


1. Introduction 


Instruments in orbit around Earth encounter a hostile evi- 
ronment of high energy photons and cosmic ray particles. 
In a Highly Elliptical Orbit a spacecraft is subjected to 
high fluxes of cosmic ray particles and gamma ray photons, 
which remain constant for long portions of the orbital 
period. For Low Earth Orbits (LEO) the natural shield¬ 
ing from charged particles due to the local geo-magnetic 


field {rigidity, see section 4.2) is much higher, but back¬ 
ground components are also highly variable on the orbital 
timescale which may be as short as an hour. For appli¬ 
cations in commercial satellites the deleterious effects of 
the background flux on electronics can be reduced by ef¬ 
fective shielding. For scientific instruments however the 
background flux is often the same energy as the photons of 
interest. In order to make calibrated determinations of the 
flux and spectra of astrophysical sources it is important 
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to accurately determine the background level. This is es¬ 
pecially true when observing the temporal characteristics 
of sources on timescales similar to that of the background 
variations. 


2. The Burst and Transient Source Experiment 


The Burst and Transient Source Experiment (BATSE) 
was an instrument aboard the Compton Gamma Ray 
Observatory (CGRO) designed primarily to detect and 
locate Gamma Ray Bursts (GRB). It consisted of eight 
modules each situated on a corner of the CGRO cuboid, as 
shown in Fig. ^ The primary instruments on each module 
were uncollimated 2025 cm^ Nal(Tl) Large Area Detectors 
(LADs), all eight of which together had a total field of 
view of 47r steradians and were sensitive from 20 keV - 2 


MeV (Fishman et al. 1989a). Data were recorded in several 
formats; the Discriminator Large Area (DISCLA) format 
contained 4 broad energy channels with data recorded ev¬ 
ery 1.024 seconds; the Continuous (CONT) data contained 
16 energy channels with a 2.048 second resolution. Each 
BATSE module also included a Spectroscopic Detector 
(SPEC) situated just below each LAD and capable of mea¬ 
suring high resolution gamma ray spectra. Lying immedi- 
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Fig. 1. The Compton Gamma Ray Observatory and one of the 8 BATSE detector modules (drawings courtesy of 
National Space Science and Technology Centre, Huntsville, AL, USA). 


ately in front of and the same size as each LAD were the 
Charged Particle Detectors (CPD), a layer of plastic scin¬ 
tillator with the main purpose of acting as a veto shield 
for the LADs. 

BATSE monitored the whole sky continuously for the 
~ 9 year lifetime of CGRO and was an important experi¬ 
ment in identifying the extragalactic origin of the Gamma 
Ray Bursts (CRB) ( Meegan et al. 1992 ). 

The 90 minute LEO of CGRO imposes a variation on 
the background which fluctuates as the spacecraft moves 
around the Earth. The background is a complicated func¬ 
tion of the pointing of the detectors, the position of the 
spacecraft in the Earth’s magnetic field and its irradiation 
history. With BATSE it is vital to know the background 
and its fluctuations in order to make accurate measure¬ 
ments of gamma ray source fluxes, when using an ob¬ 
serving technique such as the Earth Occultation Method 


Rubin et al. 

( 

1996), 

Ling et al 

( 


constructed including details of the instrument’s mechan¬ 
ical structures and the chemical composition of its compo¬ 
nent parts. A Monte Carlo numerical simulation can then 
be used to estimate the background produced when high 
energy particles interact in the detector model. A detailed 
review of mass modelling techniques, including applica¬ 
tions for BATSE and other instruments can be found in 
(| 2002 |) . 


Dean et al 


(Harmon et al. 2002). Semi-empirical models are possible, 
but they require assumptions about the physical processes 
involved or calibration data taken from a detector after 
launch and may be subject to large systematic errors. A 
different approach is used here, based on the mass mod¬ 
elling technique. 


2.1. Previous BATSE Background Models 

The estimation of the background signal in gamma ray 
telescopes has historically been a question of informed 
guesswork. Previous attempts to model the background 
count rate for BATSE detectors have been based largely 
on theoretical estimations or multi parameter fits to data 

3)). This ap¬ 


proach can now be changed due to the reduced cost of 
aquiring high performance computers and the availability 
of sophisticated particle physics code, such as the GEANT 


package (Ravndal et al. 1993). In the mass modelling ap¬ 
proach a geometrical computer model of a spacecraft is 


3. BATSE Static Mass Model 

The GEANT3 Monte Carlo particle interaction simula¬ 
tion software package has been used for several years in 
Southampton and in other groups to model astronomi¬ 
cal gamma ray detectors, suc h as the INTEG RAL Mass 
Model (TIMM), described in Lei et al. ( 1999 ). GEANT3 
is a suite of software created originally for high energy 
particle physics experiments at CERN, in which a nu¬ 
merical model of a det ector can be built an d subjected 
to fluxes of particles ( Ravndal et al. 1993 ). The suite 
of software includes detailed physical routines that use 
Monte Carlo techniques to simulate the trajectories, in¬ 
teractions and decays of sub-atomic particles within the 
model. Components within the detector are modelled as 
geometrical volumes with detailed information about their 
chemical composition (to allow accurate interaction cross 
sections and isotopic decays to be modelled by the suite of 
physics routines). The availability of computing resources 
and time limit the complexity of the model that can be 
created. Compromises can be made, however, since the 
components of most interest are those active parts of a de¬ 
tector or telescope that would record information about 
the photons or particles in the actual instrument. The 
probability of particle interaction depends approximately 
on the mass of the target volume and with each subsequent 
interaction or decay the spatial distribution of the resul¬ 
tant particles tends to isotropise. Hence the level of detail 
required for components within the model scales roughly 
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Fig. 2. BAMM contains detailed information about the BATSE modules (inset cut away image). Also included in the 
model are the structure and chemical composition of the CGRO spacecraft and the major components of the OSSE, 
COMPTEL and EGRET instruments. The lower image is a cut away along the long axis of CGRO and, to further 
optimise the view, some of the external volumes and surfaces of the model have been made transparent. 


as the ratio of the component’s mass to the square of the 
distance of that component from an active detector ele¬ 
ment. 


The BATSE Mass Model (BAMM) has been con¬ 
structed using GEANT to examine the effect of the 
space radiation environment on the detectors and has 
been adapted from earlier wo rk within the group at 
Southampton ( Dean et al. 2000 ). BAMM, shown in Fig. ||, 
is constructed from ^ 2000 geometrical volumes repre¬ 
senting the components in the entire CGRO spacecraft. 
It was possible to obtain detailed information about the 
construction of BATSE, GGRO and the other instruments 
from engineering drawings. Each BATSE detector module 
consists of ^ 100 individual volumes whose size, shape 
and chemical composition were recreated as accurately as 
possible. 


primary interest for this work are the Na(Tl) scintillator 
crystals within the BATSE LADs, although particles in¬ 
teracting in the smaller SPEG detectors and the GPDs are 
also recorded. The information recorded for any interac¬ 
tion in an active volume contains a history, including any 
previous scatterings, hadronic interactions etc., all the way 
back to the initial particle. This information is used to de¬ 
termine the initial particle’s trajectory and energetics for 
use later on in the background simulation process. Three 
main components are considered in BAMM; the cosmic 
diffuse gamma rays (CD), prompt interactions of the cos¬ 
mic rays with the spacecraft (CR) and the atmospheric 
albedo gamma rays (ATM). 

3.1. Cosmic Diffuse Gamma Ray Spectrum 


Initially the model is assumed to be in deep space, 
away from any other object, and immersed in simulated 
isotropic particle fluxes. This first stage of the modelling 
process records energy deposits from particles interacting 
in various volumes marked as ‘active’ in the model. Of 


Since there is no way of distinguishing a Cosmic Diffuse 
background gamma ray photon (CD) from a genuine 
source photon, measuring the CD background flux com¬ 
ponent is imperative in understanding the response of 
a detector. The isotropic CD background is given in 
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Fig. 3. The mechanisms behind the production of the 
atmospheric albedo gamma ray flux. Charged particles 
are indicated by solid lines whereas neutral particles are 
shown using dashes (from Zombeck (1990)). 


photons cm ^s ^MeV ^sr ^ 

(iGehrels 

1991); 

dIcD _ 
dE 

r 0.54 E-1'40 for 

0.0117 E-2-38 for 
[ 0.014 E-2-30 for 


BAMM is subjected to an 
gamma rays drawn randomly 


by the following spectra 


10 < E-y < 20 keV 

20 < < 100 keV (1) 

0.1 < < 10 MeV 

isotropic flux of simulated 
Tom the above spectrum. 


3.2. Cosmic Ray Spectrum 


The Cosmic Ray induced background component (CR) 
is simulated from an incident flux of cosmic ray protons 
given by the following spectra in protons cm“^s“^MeV“^ 
(Webber & Lezniak 1974); 


dIcR 

dEj, 


3.44 X 102f;-2 O4 

2.44 X 10 ^ 


for 3.44 <Ep< 25.5 GeV , . 
for Ep > 25.5 GeV 


The CR spectrum is assumed to be isotropic and is 
valid throughout the solar cycle in the energy ranges of 
interest. The CR component is the result of prompt inter- 
arctions in BAMM and includes gamma ray interaction 
products such as Compton scattered photons. 


3.3. Atmospheric Albedo Gamma Rays 

The Atmospheric Albedo flux of gamma rays is pro¬ 
duced by cosmic ray interactions in Earth’s atmosphere. 
A schematic diagram of the process is shown in Fig. 
Pionisation of atmospheric nuclei by the incident cosmic 
ray flux leads to the production of muons, nuclear frag¬ 
ments, and other hadrons. Gamma rays, above ~ 50 MeV 
are produced directly by the decay of neutral tt° parti¬ 
cles (rest mass ^ 140 MeV). The charged pion compo¬ 
nent also produces electrons, via the decay of charged 
muons, which radiate photons by bremsstrahlung. The 
bulk of the gamma ray flux below 50 MeV is due to the 
bremsstrahlung photons and Compton scattering of higher 



Fig. 4. A schematic illustrating the terms used in the at¬ 
mospheric albedo gamma ray emissivity function given in 
equation |[ 


energy photons. Since the photons are produced by sec¬ 
ondary processes the flux emitted by a unit volume of 
atmosphere is spatially isotropic. The energy spectrum of 
the atmospheric albedo gamma ray flux has been mea¬ 
sured accurately, with minimum contamination from the 
cosmic diffuse gamma ray background, in balloon flights 
over Palestine, Texas, USA at an atmospheric depth of 
2.5 g cm^ (pchoenfelder et a~ (1977), pchoenfelder et al 


( 1980D ). The depth dependent gamma ray flux, measured 
in photons cm“^s“^keV“^g“^ at an atmospheric depth of 


X g cm" is; 


d^ATM 

dE.y 


= 0.535xA, 


- 1.65 


( 3 ) 


Unlike the isotropic distribution of particles used when 
BAMM simulates the CD and CR components the ATM 
gamma rays are, in effect, emitted by the Earth’s disk. 
The total gamma ray emission is not isotropic over all 
zenith angles. The zenith angle distribution (measured by 
^choenfelder et ah ( 1977 )) is shown in Fig. |[ The peak at 
~ 120° is due to photons that were emitted in a preferen¬ 
tially horizontal direction (90°) but experienced Compton 
scattering between atmospheric element and detector. 


The integrated emissivity model of Dean et al. (1989) 
is used here to calculate the integrated atmospheric emis¬ 
sion of the albedo flux over the volume of atmosphere 
visible to the spacecraft, V. The line of sight dependent 
absorption between the space craft and an atmospheric 
element for photons of energy E^, characterised by the 
photon attenuation coefficient p{E~f), is given by 


M 


_ [-^(^ 7 ) J^p(r)drj 


( 4 ) 
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Fig. 5. The zenith angle distribution of the 1.5 - 10 MeV 
atmospheric albedo gamma ray flux at an atmospheric 
depth of 2.5 g cm“^ (42 km a ltitude) and a local r igidit y 
cutoff of 4.7 GV, taken from Schoenfelder et al. ( 1977 ). 
An angle of 0° in the diagram indicates downward mov¬ 
ing photons (towards the Earth) where as upward moving 
photons have an angle of 180°. The flux measurements of 
[Schoenfelder et ah (|1977|) are compared w ith those of Ling 
( 1975| ) and [Graser fc Schoenfeldei ( 1977 ). 


Hence the atmospheric gamma ray flux incident on a 
satellite is given by, 


d^ATM 

dEry 


2,h,inep(h)4E,.h]l,^cm4, 


The atmospheric air density, p, is described by two 
functions in r and h, where Fig. ^ shows the definitions of 
h, r, 9 and (/). The depth dependent mass emissivity func¬ 
tion, €[Ej, h] photons s“^keV“^g“^, is based on the spec¬ 
trum of equation ^ and takes the form; 

e(E^,x) = 0.75x°'"'E-(1-65+2-56io-^x°"“). (g) 


The emissivity function depends on the incident cosmic 
ray flux, which is shielded at low energies by the Earth’s 
magnetic held. The measurements given above were made 
over Palestine at a local rigidity of 4.7 GV. Treatement of 
the ATM component as the local value of the minimum 
geo-magnetic rigi dity c ut off changes in orbit is discussed 
further in section 4.2.2| . 

Ghanges in altitude, h, affect the solid angle subtended 
by the atmospheric disk and hence the strength of the 
albedo flux. The orientation of the detector with respect 
to the geo-centre also reduces the photon flux as the at¬ 
mosphere leaves the field of view. It would not be com¬ 
putationally possible to model every possible combination 
of altitude and pointing direction for the whole CGRO 
model. Instead a data base of responses from a single 
BATSE module is created at discrete altitudes (every 50 
km from 350 - 850 km) and attitudes (every 5° from 0 
to 180°). In the CD and CR component modelling we 
are interested not only in the direct effects of background 
photons interacting in the detector crystals, but also in 


the secondary radiation produced by high energy inter¬ 
actions in other parts of CGRO finding its way to the 
detector. However, the photons of interest in the ATM 
component are largely detected by direct interaction in 
the detectors themselves, the chance of a 2 MeV photon 
interacting in CGRO and producing a detectable photon 
in BATSE above 20 keV is small. Modelling the BATSE 
detector alone has the obvious advantage of a large reduc¬ 
tion in computing time, since the model is being run a 
total of over 400 times. 

4. Spatially and Temporally Dynamic Mass Model 

The HAMM results are ‘static’, they contain a data base 
of gamma rays and cosmic ray particles arriving from all 
directions and a data base of atmospheric albedo gamma 
rays for a grid of LAD pointing angles and orbital alti¬ 
tude. The second stage of the process is to remove those 
particles from the data base that are unphysical for partic¬ 
ular CGRO orientations and positions. The major effect 
on the background is the motion of CGRO around Earth 
modulating the component flux. The CR and CD compo¬ 
nents reduce as the Earth enters the LAD held of view 
and blocks the incident flux. The ATM component can be 
thought of as the emission of gamma rays from a spheri¬ 
cal shell and so the component flux increases as the atmo¬ 
spheric column density increases, i.e. as the limbs of the 
Earth enter the field of view. The second major effect is 
that of the local value of the rigidity of the Earth’s mag¬ 
netic field, which shields charged particles from CGRO 
and hence affects the CR and ATM components. Although 
the ATM component is effectively dependent on the CR 
flux, it is highly unlikely that a cosmic ray that interacts 
in the atmosphere and produces albedo photons will also 
go on to reach CGRO. Hence the accuracy of the model is 
not reduced by the decoupling of these two components. 

4.1. Attenuation of Components due to Obscuration 
by the Earth 

The CD and CR results from the static mass models con¬ 
tain a set of energy deposits in the BATSE LADs and 
records the momentum of the initial particle, which was 
drawn from a spatially isotropic distribution. The low 
Earth Orbit of CGRO was circular with an altitude of 
400 - 500 km. Seen from this distance the Earth’s disk sub¬ 
tended an angular radius of ^ 70° at the spacecraft, hence 
^ 30% of the sky was blocked by the Earth at any one 
time. Consequently many of the energy deposits recorded 
in the static mass model stage come from particles arriv¬ 
ing from non-physical directions and must be Altered from 
the data set. 

Since one of CGRO’s tasks was to make measurements 
of pulsars, where precise positioning is vital for barycen- 
tering pulse arrival times, the location of the spacecraft 
is accurately recorded in the telemetry data. This, along 
with the orientation of the spacecraft, allows the position 
of the Earth as seen from each of the eight LADs to be 
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calculated. Hence, given a particular part of the CGRO 
mission lifetime with its position information, the CD and 
CR components can be calculated by finding the Earth’s 
position and removing events due to Earth blockage. This 
leaves behind background components modulated for the 
effect of the Earth obscuring events from reaching the 
LADs. 


4.2. Effects of the Geo-Magnetic Field on Background 
Components 

Both the CR and the ATM components are derived from 
the incident isotropic proton cosmic ray flux. They are 
therefore both modulated, not only by obscuration by the 
Earth, but also by the shielding of charged particles by the 
Earth’s magnetic field. A cosmic ray’s ability to penetrate 
the Earth’s magnetic field is characterised by its magnetic 
rigidity, R, a measure of the particle’s momentum per unit 
charge. An ion with atomic mass number A, charge num¬ 
ber Z and kinetic energy E (in units of GeV/nucleon) has 
rigidity 

+ 2moE), (7) 


where R is measured in GV and rrio = 0.9315 GeV/c^ 
is the atomic mass unit. A particle with R less than the 
critical value, Rc, required by the local strength of the 
magnetic field will be deflected, it does not lose energy 
since the magnetic field acts in a direction perpendicular 
to the particle’s momentum. This effect shields the Earth 
from cosmic ray particles and, very simply, cuts out those 
particles with low rigidity. The rigidity required by parti¬ 
cles to reach CGRO in orbit varies as the spacecraft moves 
through the geo-magnetic field and also depends on the or¬ 
bital altitude. A simple equation can be used to estimate 
the local rigidity during a spacecraft’s orbit based on the 
Mclllwain L parameter, Rc = 15.6/ (Gehrek 1991). The 
L parameter is a coordinate describing the effective orbital 
position, as a fraction of the Earth’s equatorial radius, 
in an ideal dipole field that has the same local magnetic 
field strength as the measurement. For the circular LEO 
of CGRO at an altitude of ~ 500 km, L varies from ~ 1 
above the equator to ~ 1.5 at the highest orbital latitude 
of ±28°. This gives the minimum rigidity required by a 
proton as 12 > i?c > 7 GV, corresponding to an energy 
range of 13 < A < 8 GeV/nucleon. 

Since the GR component is modelled directly from the 
initial cosmic ray input spectrum, the BAMM results can 
be filtered according to the local minimum rigidity. The 
ATM gamma ray component is not modelled directly from 
the incident cosmic rays and so a semi-empirical approach 
must be applied. 


4.2.1. Rigidity Effects on the Cosmic Ray Component 

The shielding effect of the magnetic field can be approxi¬ 
mated by an abrupt cut-off in the incident proton cosmic 


ray spectrum, i.e. the number of particles, Nb, that would 
be detected in a local field with strength B is given by. 


Nb = 


IcndE, 


lEc 


( 8 ) 


where Ec is the minimum kinetic energy required by a 
charged particle to penetrate the magnetic field and Icr 
is the incident cosmic ray spectrum. 

In order to apply a magnetic shielding algorithm, the 
minimum local rigidity must be known. In this case a semi- 
empirical measure is made. The BATSE detector modules 
were manufactured with plastic scintillator charged par¬ 
ticle detectors (GPD) in front of each Nal LAD to act 
as a veto counter (recall Figs. and ||). The CPD rates 
were recorded by BATSE on the same 2.048 second time 
resolution as the LAD counts, positional information etc. 
Likewise, BAMM contains models of the CPD detectors 
and records particles that lead to an energy deposit within 
them in the same way that they are for the LADs. The 
BAMM CPD counts are filtered to remove those events 
that would have been prevented from reaching the instru¬ 
ment by the Earth as described in section 4.1. Assuming 
that BAMM is a reasonable estimate of the CPD angu¬ 
lar response, the number of charged particles remaining is 
that that would be detected by the CPDs in the absence of 
any magnetic field. Using equation ^ an approximation for 
the ratio of CPD counts observed by BATSE, VcpDobs, to 
those predicted by the model, VcpDmod, can be made; 


A^CPPobs _ Iec ^cndE 
A^CPDmod /q°° IcndE 

Solving the integral formed by substituting the cosmic 
ray spectrum and boundaries given by equation || into 
equation gives an expression for VcpDobs/A^CPDmod in 
terms of Ec. This can be rearranged to give an expression 
for Ec] 


Ec = 


0.0601 




V NcPDmod 


- 3.743 X 10 


-3 


329.89 


( 10 ) 


Equation can be used to determine the value of Rc 
from Ec. At each 2.048s time interval, for which the po¬ 
sition information for CGRO is known, the value of Ec is 
calculated and compared against the incident particle en¬ 
ergy for each energy deposit predicted by BAMM. If the 
incident particle is determined to have energy less than 
this minimum value it can never have reached a position 
to cause an energy deposition and it is removed from the 
data base. 


4.2.2. Rigidity Effects on the Atmospheric Component 

Background components that are dependent on the cosmic 
ray flux but are not modelled from individual incident 
cosmic rays, such as the ATM gamma ray component, 
must be treated differently. 
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An empirical rigidity correction factor, Af(i?c), based 


on a number of satellite observations, e.g. Imhof et al. 
(1976), may be defined as. 


M{R) = 



dljEp) 

dEp 


E°-'^dE^ 


( 11 ) 


where dI{Ep)/dEp is the cosmic ray intensity at pro¬ 
ton energy Ep. Numerous measurements have been made 
of Rc above Palestine, Texas, of 4.7 GV. A correction fac¬ 
tor, f{R) = M{Rc)/M{Rc = 4.7), is defined that can be 
used to scale the Palestine rigidity measurements to any 
value of Rc- By solving the integral formed by substituting 
equation ^ with the cosmic ray spectrum and boundaries 
of equation the rigidity correction factor is, 

f{R) = 28.462A-°'3^2 _ Q 579^ (72) 


channels by individual detectors during periods of the 
orbit when they are pointing towards the poles. The ef¬ 
fect manifests itself as a short ‘blip’ - a jump of a few 
tens of counts that lasts ^ 5 minutes and appears two or 
three times a day in detectors alligned in favourable direc¬ 
tions. Since these events are short and infrequent they are 
not considered a serious problem. Modelling the electron 
bremsstrahlung component would also be difficult since 
the strength of the effect is dependent on solar activity 
amongst other things. 

A detailed study of these two effects has shown that 
a loss of only ^ 15% of usable data are achieved by not 
considering parts of the dat a that are s e verely affecte d 
by these two components (Throp 2001; Brown] ^OOlj ). 
Although it is in principle possible to use the mass mod¬ 
elling technique to produce these extra components it is 
not considered an efficient use of time for the present. 


4.3. Other Components 


4.4. Renormalization of Background Components 


The background components mentioned above are mod¬ 
ulated on a psuedo periodic time scale, governed by the 
orbital motion of CGRO around the Earth, and are rea¬ 
sonably simple to predict. Other particle processes ex¬ 
ist that can produce a detectable background component 
with weak dependence on the orbital period. These back¬ 
ground contributions may become important for specihc 
spacecraft positions and/or orientations and as such are 
harder to predict than the CD, CR and ATM components. 
Two effects, due to the South Atlantic Anomaly (SAA) 
and atmospheric electron bremmstrahlung are considered 
here: 

The very high fluxes of energetic protons found within 
the SAA (> 1000 protons cm“^s“^ above 30 MeV) acti¬ 
vate isotopes within the material of a sp acecraft passing 
through it (see for example Roederei 1970| ). These isotopes 
subsequently decay and produce a delayed count rate in 
gamma ray detectors that can last for tens of minutes af¬ 
ter the initial pass. The charged particles are also capable 
of severely damaging electrical equipment passing through 
them. A typical solution, used with BATSE and CGRO, is 
to avoid damage by shutting down instruments during pas¬ 
sage through the SAA (about 5-6 times during a typical 
24 hour period). This makes the estimation of the decay of 
the spectral count rates on emergence from the SAA un¬ 
reliable. The problem is further compounded since much 
of the incident particle flux has low penetrating power in 
the material of the spacecraft. Therefore anisotropy of the 
trapped particle momentum distribution with respect to 
the CGRO orientation vector makes exact determination 
of the induced radioactivity, and hence the eventual back¬ 
ground contribution, extremely difficult to model. 

The Earth’s magnetic field is also responsible for 
a small background component caused by aurora¬ 
like emissions of x-ray photons produced by elec¬ 
tron bremsstrahlung above the North and South poles 
(Aschwanden 1998| ; Aschwanden et al. 1998] ). Eor BATSE 
these emissions may be observed in the lowest energy 


The mass modelling technique outlined above produces 
three components, which when added together, gives a 
close approximation of the actual gamma ray background. 
The accuracy of the technique is limited by how well the 
input spectra for the CD, CR and ATM components have 
been measured. Due to uncertainty in the shape of the 
CD, CR and ATM input spectra the results of the mass 
modelling cannot exactly reproduce the actual count rates 
seen in the BATSE detectors. Although the time profile 
of the counts seen in each detector is reproduced accu¬ 
rately the absolute normalization of each component is 
not. Consequently it is necessary to rescale the three com¬ 
ponents. 

A short section (^ 5 orbits) for each detector and 
CONT energy channel in a time period of interest is 
chosen such that it is free from contamination due to 
SAA passages, known flaring gamma ray sources and large 
data gaps. This section is then corrected by subtracting 
a contribution from all known gamma ray sources above 
200 mCrab brightness. A model of the form 

C(t) = a • CD(t) -k 3 ■ CR(t) + 7 • ATM(t) (13) 

is fitted to the short section of corrected data using 
a least squares method. The rescaling returns the global 
amplitudes, a, /3, 7 , of the CD, CR and ATM components 
respectively which are used to return the total background 
for the day, detector and channel in question. It is impor¬ 
tant to highlight that the htting routine in no way alters 
the predicted shape of the count rate profile of the three 
components returned by the mass model. 

5. BATSE Dynamic Mass Model Results and 
Applications 

Figure ^ shows the excellent agreement between the dy¬ 
namic BAMM data and the CONT count rate profile in 
the BATSE LADs for three energy channels. The data 
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Fig. 6. BATSE LAD 5 dynamic background model for three different energy channels in a twelve hour section of the 
CGRO mission. The top row of figures shows the excellent agreement of the dynamically filtered BAMM (dark line) 
with the BATSE CONT data (light dots). The central figures show the residuals of the above plots and highlight the 
large reduction in the BATSE background fluctuation that can be achieved by subtracting the model. The bottom 
diagrams show the CD (black), ATM (mid grey) and CR (light grey) components of the total background model and 
a DC offset due to unresolved background gamma ray sources (dotted line). The dependence of the three components 
on energy can clearly be seen. 


are flat-fielded simply by subtracting the modelled back¬ 
ground, which reduces the fluctuation of the background 
counts by up to an order of magnitude. The relative im¬ 
portance of the three components as a function of energy is 
clear; at low energies the CD component dominates the to¬ 
tal background; at the high energy end the CR component 
is most important; the ATM component is at its strongest 
at intermediate energies. The figure also shows the impor¬ 
tance of modelling the detector attitude correctly. The CD 
component shows a large modulation, indicating that the 
detector was pointing almost directly in the plane of the 
CGRO orbit. As the pointing angle moves away from the 
orbital plane less of the Earth comes into the detector field 
of view and so less of the component will be blocked out. 
In this case the component will have a similar maximum 


strength but much reduced variation. Two applications of 
the flat-fielded data are discussed below. 


5.1. Gamma Ray All Sky Survey 


The major driving force behind developing the dynamic 
mass model for BATSE was to enable the systematic er¬ 
rors in long (~ hundreds of days) sections to be reduced 
such that it could then be used for continual observa¬ 
tions of objects with improved sensitivity. The most im¬ 
portant application of this is in the production of a hard 
x-ray all sky survey using the Earth occultation technique 
(Fishman et al. 1989b| ; Barmon et al. 2002). A maximum 
likelihood method is used to produce all sky images in the 
gamma ray band, the first time this has been achieved 
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Fig. 7. Maximum likelihood all sky image, displayed in 
Galactic Coordinates and showing the effect of applying 
the LIMBO technique to raw data. The image is made 
from CONT energy channel 2 (~ 35 - 45 keV) for the pe¬ 
riod TJD 11400 - 11450 and the peak flux significance in¬ 
dicated is at a position consistent with the Crab Nebula, 
(/,6) = (-175.4,-5.8). Note the large systematic back¬ 
ground, present at a level corresponding to ^ 25% of the 
Crab flux. 


since the HEAO missions of the late seventies, described 


Levine et al. (1984). The Likelihood Imaging Method 


for B ATSE Occul t ation data (LIMBO) are discussed in 
full in Shaw et al. ( 2001 , ^002 ). 



35 keV - 45 keV Flux Significance limited to max of 0.50 * 89.27 a 
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Fig. 8. Maximum likelihood all sky image, as Fig. |^, but 
applied to data flat-fielded with dynamic BAMM data. 
The same greyscale has been used in both figures. The 
systematic background apparent in Fig. ^ has been re¬ 
duced significantly. 


The flat-fielded data was shown to have much re¬ 
duced background variation in figure and thus should 
produce much cleaner images. Shown in Figs. and 
are the results of applying LIMBO to the raw and flat- 
fielded BATSE count rates respectively. The maps were 
created for one CONT energy channel (~ 35 - 45 keV) and 
summed for the 50 day period from TJD 11400 - 11450. 

A severe and varying systematic background is present 
across the map made from the raw data. Near the Crab 
Nebula position it is at a level equivalent to ~ 30(T or 
^ 250 mCrab. In contrast the flat-fielded map shown in 
Fig. H has a drastically reduced systematic level that, at 
worst, is ~ 50 mCrab. By simple extrapolation of the sig¬ 
nificance of the measurement of the Crab Nebula flux (but 
including any systematic offset) it is possible to estimate 
the statistical 3 a flux limits for the maps: ~ 25 mCrab 
and ~ 30 mCrab respectively for the raw and flat-fielded 
images. The systematic error is clearly the dominant fac¬ 
tor in determining the sensitivity of the imaging process. 
When the entire BATSE data base is considered we ex¬ 
pect to produce an all sky survey map with a sensitivity 
of ^ 1-2 mCrab. 

The faint features stemming from bright point sources 
in the flat-fielded map are due to the imaging process it¬ 
self. Individual point sources are created by the superpo¬ 
sition of many images of the Earth’s limb (radius ^ 70°). 
Therefore it is possible for a bright source to have an ef¬ 
fect many tens of degrees away from its position on the 
sky. The ‘point spread function’ of a single source is a 
complicated function of many parameters, but it is calcu¬ 
lable. Further reduction of the systematic background in 
the images, by producing methods of image cleaning the 
occultation maps and by improving the background mod- 
elling , remains ongoing work ( Shaw et al. 2002 ; Westmore 

S))- 


5.2. Gamma Ray Burst Studies 

In BAMM the accurate geometrical and chemical rep¬ 
resentation of the BATSE modules and their immedi¬ 
ate neighbourhood (i.e. the entire CGRO spacecraft) may 
have further useful applications in the analysis of Gamma 
Ray Bursts (CRB). The GRBs detected by BATSE emit 
up to 10^'^ ergs above ~ 30 keV (postulated for isotropic 
emission) over a range of times as short as a few hun¬ 
dredths up to ^ 1000 seconds ([Fishman fc Meeg^ 1995 


Meszaros 1997). Despite the higher signal to background 
ratio seen in CRB detections than in other types of gamma 
ray sources a background model is still important since the 
fluxes in CRB lightcurves can vary greatly. It is also tan¬ 
talising to be able to investigate the most distant GRBs 
that are detected at the threshold of sensitivity. 

For the most part the short length of the burst is such 
that the background does not vary appreciably and can 
be estimated by fitting simple empirical models to the 
data (e.g. Preece et al. ( pOOO )). However for the fraction 
of GRBs that last for more than a few 100 seconds these 
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Fig. 9. The modelled response of one of the BATSE SPEC 
spectrometer detectors to a GRB. In this simulation the 
input spectrum was based on the spect rum m e asured 
for GRB91053 (Parameters provided by Briggg (^001 ) 
and using the model of [Band et al. ( 1991 ): A = 0.118 
photons cm“^s“^keV“^, Eo = 753 keV, a = —0.907, 
(3 = —2.72).The figure suggests that 10-30% of the to¬ 
tal detected flux is not detected directly from the gamma 
ray burst but is in fact due to secondary photons emitted 
when gamma ray photons interact in the material of the 
space craft. 


models start to run into difficulties as the background 
starts to vary. In addition gaps in the BATSE telemetry 
become more likely the longer the time being considered 
and this further compounds problems with the fit. 

On timescales of many hours to perhaps a few weeks 
the meas urement of GRB tails has become one of great 
interest ( Giblin et al, 1999| ). The current method of esti¬ 
mating the background for these events has involved inter¬ 
polation of BATSE data from orbits before and after the 


orbit of interest (Connaughton 2002). This method has 


shown interesting results, but is subject to systematic er¬ 
ror in determining the count rate due to the precession of 
CGROs orbit and can be affected by change in the space¬ 
craft conditions such as a reorientation to a new pointing. 

Spectral studies of GRBs is also an area of high sci¬ 
entific interest and a great deal of effort has been put in 
to searching for lines and other features in BATSE GRB 
spectra (e.g Briggs et al. ( 1998 ), Palmer et al. ( |1994 ) and 
others). Knowledge of the spectral response of the detec¬ 
tors is highly important in deconvolving GRB spectra. 
However, this is difficult to estimate due to the partial 
deposition of energy by high energy photons in the de¬ 
tector material (Briggs 1996) and the detection of sec¬ 
ondary radiation created by Compton scattering of the 


source photons in the neighbourhood of the detectors, such 
as the spacecraft material or Earth’s atmosphere (Hua 
fc Lingenfelter 1996|). Pendleton et al. (19951) have used 


a Monte-Carlo model to estimate the Detector Response 
Matrices (DRM) of the LAD and SPEC detectors. This 
model assumed azimuthal symmetry of the detector re¬ 
sponse, which has limited the DRMs use with SPEC data 
to events occuring at no more than 60° from the detector 
axis. The representation of the surrounding CGRO instru¬ 
ment was also fairly crude and the model is not able to 
account for the effect in all eight of the detector modules 
together and in situ on CGRO for different orientations 
and in different magnetic field conditions. 

BAMM may be a very useful tool in future GRB analy¬ 
sis. With a rigourous treatment of the effects on the spec¬ 
trum due to the surrounding material, it will allow the 
measurement of the gamma ray background in all eight 
BATSE detectors for any position and orientation in the 
orbit history on any timescale down to the maximum time 
resolution of the dataset. As an example Fig. ^ shows the 
spectrum calculated from one of the SPEC detectors by 
BAMM when presented with a simulated flux of gamma 
ray photons from a GRB. A significant fraction (^ 10 
- 30 %) of the total flux is due to secondary photons 
and would lead to serious systematic errors when calcu¬ 
lating the spectrum of the GRB. The high level of spa¬ 
tial, mechanical and chemical information within BAMM 
together with the particle interaction simulation capabil¬ 
ities of GEANT3 also allows particles to be tracked in 
detail through the detector. For example it is possible for 
BAMM to identify the high energy particles that deposit 
a small amount of energy, leave the detector and are then 
detected again after further scattering interactions else¬ 
where in the instrument. Work continues to investigate to 
what level BAMM can be used and how the mass mod¬ 
elling technique can be adapted to other missions, such 
as the BAT GRB monitor onboard the SWIFT m ission 
(Willis 2002; Caraveo fc Swift Collaboratioi] 2001). 
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